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ABSTRACT 

We present a set of Roche tomography reconstructions of the secondary stars in the 
cataclysmic variables AM Her, QQ Vul, IP Peg and HU Aqr. The image reconstruc- 
tions show distinct asymmetries in the irradiation pattern for all four systems which 
can be attributed to shielding of the secondary star by the accretion stream/column in 
AM Her, QQ Vul and HU Aqr, and increased irradiation by the bright spot in IP Peg. 
We use the entropy landscape technique to derive accurate system parameters (Mi, 
M2, i and 7) for the four binaries. In principle, this technique should provide the most 
reliable mass determinations available, since the intensity distribution across the sec- 
ondary star is known. We also find that the intensity distribution can systematically 
affect the value of 7 derived from circular orbit fits to radial velocity variations. 

Key words: binaries: close - novae, cataclysmic variables - stars: late-type - stars: 
imaging - line: profiles. 



1 INTRODUCTION 

The secondary, Roche-lobe-filling stars in CVs are key to 
our understanding of the origin, evolution and behaviour 
of this class of interacting binary. To best study the sec- 
ondary stars in CVs, we would ideally like direct images of 
the stellar surface. This is currently impossible, however, as 
typical CV secondary stars have radii of 400 000 km and 
distances of 200 pc, which means that to detect a feature 
covering 20 per cent of the star's surface requires a resolu- 
tion of approximately f microarcsecond, f 000 times greater 
than the diffr action-limited resoluti on of the world's largest 
telesco pes. iRutten fc Dhillonl Jf994l) and lWatson fc Dhillon] 
teOOll hereafter referred to as Paper I) described a way 
around this problem using an indirect imaging technique 
called Roche tomography which uses phase-resolved spectra 
to reconstruct the line intensity distribution on the surface 
of the secondary star. 

Obtaining surface images of the secondary star in CVs 
has far-reaching implications. For example, a knowledge of 
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the irradiation pattern on the inner hemisphere of the sec- 
ondary star in CVs is essential if one is to calculate stellar 
masses accurately enough to t est binary star evolution mod- 
els fsee lSmith fc Dhillonlf993l . Furthermore, the irradiation 
pattern provides information on the geome try of the acc ret- 
ing structures around the white dwarf (see lSmithllf 995lb 

In this paper we present new Roche tomograms of the 
secondary star in the dwarf nova IP Peg, and the mag- 
netic CVs AM Her and QQ Vul in the light of the Na 
I AA8I83,8195 A absorption doublet. In addition, we also 
present a Roche tomogram of the magnetic CV HU Aqr 
in the light of the He II A4686 A emission line component 
known to originate from the secondary star. These tomo- 
grams allow a study of the irradiation pattern on the sec- 
ondary stars as well as providing measurements of the binary 
parameters for all four CVs. 



2 OBSERVATIONS AND REDUCTION 

The spectra of QQ Vul and HU Aqr were taken 
on the 3.5-m telescope at Calar Alto and full de- 
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tails of the obser vations and subsequent da t a redu ction 
can be found in ICatalan. Sc hwope fc Smith! il999T) and 
ISchwope. ManteTA^HoTn^n^fflTf) . respectively. 

The observations of IP Peg and AM Her were carried 
out over two nights on 1994 August 11-12 using the ISIS 
dual-beam spectrograph and 1200 lines mm -1 gratings on 
the 4.2m William Herschel Telescope (WHT). The EEV3 
CCD chip in the red arm covered the wavelength range 
7985-8437A at a resolution of 0.37A (~27 km s" 1 ). The 
TEK2 CCD chip on the blue arm covered the wavelength 
range 4560-4968A at a resolution of 0.4A (~50 km s" 1 ). 
To reduce readout noise, the pixels were binned by a factor 
of 2 in the spatial direction. Each object was observed using 
exposure times of 300 seconds and a nearby field star was 
located on the slit to correct for slit losses. Wide slit spectra 
of the companion star were also taken to permit absolute 
spectrophotometry. Sky-flats and tungsten flats were taken 
on both nights and a CuAr+CuNe arc spectrum was taken 
approximately every 50 minutes in order to wavelength cal- 
ibrate the spectra. Narrow and wide slit spectra at a range 
of different air-masses were taken of flux standards for the 
purposes of flux calibration and water removal. We also ob- 
served a range of M-dwarf spectral-type templates in order 
to obtain a Na I intrinsic line profile for use in Roche to- 
mography. 

During the reduction process, bias subtraction, division 
by a master flat-field frame and sky subtraction were carried 
out a nd then the spectra were optimally extracted llHornel 
1986). The arc spectra were used to wavelength calibrate 
the data, with an rms error of O.OlA. Slit loss correction 
was then applied by dividing the object spectra by the av- 
erage value per pixel (after masking out strong lines, cosmic 
rays and telluric features) in the corresponding comparison 
star spectra, and later multiplying by the average value per 
pixel in the wide-slit comparison star spectrum. Telluric ab- 
sorption and absolute flux calibration were then performed 
using the standard star spectra. 

In total we obtained 79 spectra of IP Peg and 51 spec- 
tra of AM Her. We present only the red spectra in this 
paper, as we are only concerned with the Na I doublet at 
AA8183,8195A for the purposes of this work. A full journal 
of observations is presented in Tabled 



3 ROCHE TOMOGRAPHY 
3.1 Techniques 

Roche tomography dRutten fc Dhillonl Il994t 

iDhillon fc Watsonl yOOl; Paper I) is analogous to the 
Dopplcr imagin g technique used to m ap rapidly rotating 
single stars (e.g. IVoet fc Penrodlll983l) . In Roche tomogra- 
phy the secondary star is assumed to be Roche-lobe-filling, 
locked in synchronous rotation and to have a circularized 
orbit. The secondary star is then modelled as a series of 
quadrilateral tiles or surface elements, each of which is 
assigned a copy of the local (intrinsic) specific intensity 
profile. These profiles are then scaled to take into account 
the projected area of the surface element, limb darkening 
and obscuration, and then Doppler shifted according to 
the radial velocity of the surface element at that particular 
phase. Summing up the contributions from each element 



gives the rotationally broadened profile at that particular 
orbital phase. 

By iteratively varying the strengths of the profile con- 
tributed from each element, the 'inverse' of the above pro- 
cedure can be performed. Due to the variable and unknown 
contribution to the spectrum of the accretion regions in 
CVs, the data must be slit-loss corrected and continuum 
subtracted prior to mapping with Roche tomography. Thus, 
Roche tomograms present images of the distribution of line 
flux on the secondary star. The data is fit to an aim and 
a unique map is selected by employing t he maximum en- 
tropy (MEMSYS) algorithm developed bv lSkilling & Brvanl 
( 1984). The definition of image entropy, S, that we use is 

k 

S =y rrij- dj - m 3 In (rrij/dj) 
j=i 

where k is the number of tiles in the map, and rrij and dj 
are the map of line flux and the default map, respectively. 
We employ a moving uniform default map, where each ele- 
ment in the default is set to the average value in the map. 
For furth er details of the Roche tom ography process see the 
review bv lDhillon fc Watsonl (1200 ill . 

For those CVs where we wish to map the Na I doublet, 
Roche tomography requires the intrinsic line profile to be 
input as a template. Ideally, this should be the same line 
taken from the spectrum of a slowly rotating star of the 
same spectral class as the star being mapped, preferably 
obtained with the same instrumental set-up. Unfortunately, 
only one of the spectral-type template stars that were ob- 
served is useful (GL806 - see Tabled; the others show signs 
of broadening of the Na I doublet and have been identified 
as either flare stars, and are therefore probably rapidly ro- 
tating, or binary stars, where V ro t sin i measurements using 
the Na I absorption doublet have been found to be system- 
atically greater tha n those given by the narrow metal lines 
jBleach et alJl200Gft . We used the GL806 template in all of 
our reconstructions and, even though the spectral class of 
the template differs from those of the stars being mapped, 
the exact shape of the template profile used has little effect 
on the reconstructions due to the large degree of rotational 
broadening present in CV secondary stars. 

For the CVs mapped using the Na I doublet, it is nec- 
essary to correct the derived systemic velocities for the sys- 
temic velocity of the spectral-type template used in the re- 
constructions. The systemic velocity of GL806 was deter- 
mined by fitting gaussians to the Na I doublet, which yielded 
a valu e of 7 = -24.3 ± 1.1 km s _1 . For comparison. IWilsonl 
lll953l) found the systemi c veloc ity of GL806 to be -15 ± 
5 km s _1 , whereas iGliesd lll969T) quote a value of -24.4 km 

It is important to note that the maximum-entropy al- 
gorithm we use requires that the input data are positive. 
As a result, it is necessary to invert absorption-line profiles 
and, in the case of low signal-to-noise data, add a positive 
constant in order to ensure that there are no negative values 
in the continuum-subtracted data. If this positive constant 
is not added, and any negative points are either set to zero 
or ignored during the reconstruction, then the fit will be 
positively biased in the wings of the line profile wings. The 
exact value of this constant is unimportant, so long as it 
prevents the occurrence of negative values in the data and 
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Table 1. Jo urnal of observatio ns. Orbital phase is shown for AM Her using the ephemeris of lSouthwell et all and IP Pee; using the 

ephemeris of Wolf et alJ (1993J). Since the latter ephemeris is deduced from the white-dwarf egress in IP Peg, an appropriate correction 
has been made to ensure that phase corresponds to the superior conjunction of the primary. Note: HU Aqr was trailed along the slit 
during the exposures. The final time resolution was ~ 30 s, and the data was phase-folded into 100 phase bins covering a complete orbital 
cycle. 



Object 


UT Date 


UT 


UT 


Phase 


Phase 


± cxp 


No. of 




(D/M/Y) 


start 


end 


start 


end 


00 


Spectra 


AM Her 


11/08/94 


21:38 


23:58 


0.57 


1.29 


300 


22 


AM Her 


12-13/08/94 


21:20 


00:09 


2.22 


3.11 


300 


29 


QQ Vul 


08-09/07/91 


23:33 


04:10 


0.12 


1.33 


345-720 


16 


QQ Vul 


09-10/07/91 


21:47 


03:52 


6.12 


7.71 


600-680 


20 


IP Peg 


12/08/94 


01:23 


05:26 


0.88 


1.93 


300 


41 


IP Peg 


13/08/94 


01:13 


05:11 


7.16 


8.18 


300 


38 


HU Aqr 


17-18/08/93 










750-1800 


7 



is effectively removed. We account for this positive constant 
during the reconstruction by employing a virtual image ele- 
ment which contributes a single value to all data points, ef- 
fectively cancelling out the constant. For the systems where 
the Na I doublet has been used in the reconstructions, we 
have employed this technique of adding a positive constant 
to each data value. We find that this significantly improves 
the reliability of the technique to determine the component 
masses when noisy data are being used. As the data for HU 
Aqr are from Gaussian fits to the He II A4686A emission 
line, there are no errors in the continuum and hence we do 
not apply this technique. 

In addition, we have accounted for the effects of velocity 
smearing due to the motion of the secondary star over the 
duration of an exposure for all four systems. This is done by 
calculating the profile at evenly spaced intervals over one ex- 
posure and then averaging them. Although computationally 
more time-consuming, this reduces the effects of smearing, 
and can also take into account features that may cross the 
limb of the secondary star during an exposure. 

3.2 Determining the system parameters 

Roche tomograms are greatly influenced by parameters such 
as the component masses, inclination and systemic velocity 
of the binary. Errors in the assumed values of these param- 
eters degrades the quality of the surface map and, in gen- 
eral, results in additional structure in the Roche tomograms. 
The correct parameters are, therefore, those that produce 
the map containing the least artefacts, corresponding to the 
map of highest entropy (Paper I). By carrying out recon- 
structions for many pairs of component masses (iterating to 
the same % 2 on each occasion) we can construct an 'entropy 
landscape' (e.g., Fig.[3J, where each point corresponds to the 
entropy value obtained in a reconstruction for a particular 
pair of component masses. 

In order to determine the binary parameters of the CVs 
studied in this paper, we constructed a series of entropy 
landscapes assuming a wide range of different values for the 
orbital inclination and systemic velocity. For each combina- 
tion of i and 7, we picked the maximum entropy value in the 
corresponding entropy landscape, and plotted this value in 
the i — 7 plane, which is shown for each object on the left- 



Table 2. Spectral- type templates observed with the WHT. 



Object 


Spec. Type 


Notes 


GL 65B 


M5.5V 


flare star 


GL 866 


M5.5V 


flare star 


GL 83.1 A 


M4.5V 


flare star 


GL 699 


M4Ve 


variable star 


GL 725A 


M3V 


double star 


GL 806 


M1.5V 


used in this paper 



hand side of Fig. Q In all cases we found that a unique or 
optimal systemic velocity could be selected that consistently 
gave the map of maximum entropy. In addition, the value 
of the optimal systemic velocity was found to be largely 
independent of the orbital inclination assumed during the 
reconstructions, as demonstrated in the left-hand plots of 
Fig. Q For example, in reconstructions carried out over a 
range of inclinations spanning 30° , the systemic velocity we 
obtain for both IP Peg and AM Her only varies by 1 km 

The inclination, however, is not as well constrained for 
any of the CVs. The right-hand plots on Fig. Q show cuts 
through the plots in the i — 7 plane at the optimal systemic 
velocity, and the scatter clearly demonstrates the difficulty 
in selecting an optimal inclination, although the deteriora- 
tion in the quality of the reconstruction for some inclinations 
allow limits to be set. 



The difficulty in constraining the inclination can be ex- 
plained by considering that there are two ways of determin- 
ing the inclination of non-eclipsing CVs from the spectra of 
the secondary star alone. First, from the variation in the 
projected radius of the Roche-lobe shaped secondary star 
as it is viewed at different aspects, corresponding to a vari- 
atio n in the shape (or measured usini) of the line profiles 
(see Shahbaz 1998). The second is from the variation in the 
strength of the profiles, since a high inclination system will 
exhibit a much greater variation in the line profile strength 
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Table 3. Adopted system parameters. 





AM Her 


QQ Vul 


IP Peg 


HU Aqr 


i 


60-80° 


72° < i <90° 


82°- 85° 


84.4° 


1 


0.53 


0.63 


0.43 


0.25 


Mi (M ) 


~0.36 


0.58-0.66 


1.16-1.18 


0.61 


M 2 (M ) 


~0.18 


0.34-0.44 


0.50 


0.15 


7 (km s~ 1 ) 


-17 


-14 


-9 


+7 



than a low inclination system with the same surface intensity 
distribution. 

Unfortunately, the variation in v sin i is expected to be, 
at most, around 20-25 km s" 1 for high inclination systems. 
Although this is comparable to the velocity resolution of our 
data, it is buried in noise and hence the mapping technique 
is largely blind to any such variations. Our only constraint 
on the inclination, therefore, comes from the variation of 
the line strength, which also has its limitations. Although, 
rapidly varying features can only be mapped onto high in- 
clination systems, slowly varying features produced by a low 
inclination system can also be mapped onto the polar regions 
of a high inclination system with little loss in the quality of 
the reconstruction. Thus we only expect to constrain the in- 
clination in CVs which have high inclinations and exhibit 
strong asymmetries in their surface intensity distributions. 
In addition, for the same reason outlined above, we also 
expect our determinations to be biased towards higher in- 
clinations. 

As we cannot reliably constrain the inclination of the 
non-eclipsing CVs in this work, and hence the component 
masses, Fig. shows the derived component masses as a 
function of inclination. The top section of each panel in 
Fig. [5] also shows the mass ratio derived for each inclina- 
tion. Since vsini is purely a function of the mass ratio q 
(=Mi/M2) and the radial velocity semi-amplitude of the 
secondary star K r , reconstructions over a range of inclina- 
tions should yield the same value of the mass ratio on each 
occasion. This is confirmed in Fig- 01 with the scatter in the 
values serving as an indication of the noise in the technique. 
Indeed, the scatter can be accounted for by the fact that the 
component masses are only varied in O.O2M0 increments in 
the entropy landscapes. 

Entropy landscapes representative of the mass esti- 
mates for each system are shown in Fig. |3] and a summary 
of the adopted system parameters can be found in Table [3] 
The results for each CV are discussed in more detail in the 
relevant sections below. 



and -40 ±15 km s 1 from circular fits to the H7 and He I 
A447lA emission lines, respectively. 

Fig. also shows the maximum entropy value as a 
function of inclination obtained in entropy landscapes con- 
structed assuming an optimum systemic velocity of -17 km 
s _1 . It appears as though we cannot determine a unique 
inclination for AM Her, though our work suggests that 
the inclination may lie around 60-80°. Given that we ex- 
pect our results to be biased towards higher inclinations, 
the inclination is most likely to be around 60°. Although 
this estimate is far grea ter than the ~35° determined by 
Braincrd <§^LambJ (U985), it is consistent with the findings 
of lDavev fc Smith! il996T> who used radial- velocity curve fit- 
ting to constrain t he in clination to be between 45° and 
60° . iGreelev et alJ (I1999T) observed a partial eclipse of the 
He II A1640 A emission line region by the secondary star 
and hence inferred an i nclina tion >45°. In addition to this, 
IWickramasinghe et alJ lll99ll) obtained i = 52° from polari- 
metric observations. 

Fig. [3] shows the entropy landscape for AM Her con- 
structed using i — 60°, 7 = -17 km s _1 and taking into 
account velocity smearing introduced by the 300 second ex- 
posure times. The brightest point on this plot indicates the 
map of highest entropy, and corresponds to Mi = 0.36 M 
and M2 = 0.18 M , which we adopt as our best estimate 
of the component masses. The optimum component masses 
decrease to Mi = 0.24 M and M 2 = 0.12 M for i = 80°. 
These masses appear suspiciously low since AM Her lies 
above the period gap and we would not, therefore, expect the 
secondary mass to be significantly less than AI2 = ~0.25 M . 
This partly confirms our speculation that our inclination 
determinations are biased towards high inclinations (Sec- 
tion and that the inclination is not likely to be greater 
than 60°. 

The mass of the secondary star derived from the en- 
tropy landscape agrees well with the value of M2 = 0.20 
- 0.26 M estimated by ISouthwell et alJ (I1995T) . This also 
agr ees with the main - seque nce mass-period relation derived 
bv ISmith fe DhillorJ (1998), which gives a secondary star 
mass of 0.23±0.02 M . 

Meaningful comparison of the primary mass derived 
in our work with other published white dwarf masses 
is, unfortunately, more difficult as they cover a wid e 



range of v alues, including: 0.39 M 1 Yoim g et alJ ll98lTl , 



!).!,'.! M^ iWu. Chanmugam fc Shavivl 11993) . 0.75 M 



jMukai fc CharlesHl987T). 0.91 M jMouchetlll993fo arid 1.2 2 
M^ dCropper. Ramsav fc Wulll998l) . lGansicke etaH <199& . 

however, concluded that the white dwarf mass in AM Her 
was between 0.35 - 0.53 M and ruled out a white dwarf 
mass as high as 1.22 M . In addition to t heir secondary 
star mass estimate, ISouthwell et alJ il995T> derived q = 
0.47±0.08, in agreement with the value of q = 0.53 obtained 
in this work (Fig. |2J . 



4 RESULTS FOR AM HER 

4.1 The binary geometry of AM Her 

From Fig. the maximum entropy value yields a value 
of -17 km s _1 for the systemic velocity of AM Her, tak- 
ing into account the systemic velocity of the template 
star. This is within 2a of th e rang e of values found by 
lYoung. Schneider fc Shectmanl lll98ll) of -14 ± 4 km s _1 



4.2 The surface map of AM Her 

The Roche tomogram for AM Her is presented in Fig.^Jand 
shows a decrease in the amount of Na I absorption around 
the inner Lagrangian (Li) point on the trailing hemisphere. 
We tested the significance of this feature using the following 
procedure: 
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Figure 1. The plots on the left show the maximum entropy value obtained in each entropy landscape when a particular combination 
of inclination and systemic velocity is employed during the reconstructions; horizontal cuts through these plots at the optimal systemic 
velocity found for each object are shown on the right. From top to bottom: AM Her, QQ Vul, IP Peg and HU Aqr. The solid curves 
show parabolic fits to the data points. 
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Figure 2. Estimates of the component masses for different inclinations. Clockwise from top left: AM Her, QQ Vul, HU Aqr and IP Peg. 
The top panel of each plot shows the mass ratio, the lower panel shows the individual masses with the upper curve showing the primary 
mass and the lower curve showing the secondary mass. 



(i) Fit the 'true' map using the observed data (Fig.^J. 

(ii) Compute 200 'trial' maps from simulated data sets 
constructed using bootstrap re-sampling. From our observed 
trailed spectrum containing n data points we formed a sim- 
ulated trailed spectrum by selecting, at random and with 
replacement, n data values and placing these at their orig- 
inal positions in the new simulated trailed spectrum. For 
points that are not selected, the associated error bar was set 
to infinity and hence they were effectively omitted from the 
fit. For points selected more than once, the error bars were 
divided by the square root of the number of times they were 
picked (see Paper I for more details). 

(iii) Calculate 200 'difference' maps by subtracting each 
'trial' map from the 'true' map. These 'difference' maps give 
the scatter in each pixel. (It would be incorrect to now cal- 
culate a summary error statistic like the standard deviation, 
as the distribution of pixel values is often found to be non- 
normal; see Paper I). 

(iv) Subtract some mean or comparison level from the 
'true' map. For instance, if one wished to determine the sig- 
nificance of a starspot, an appropriate comparison level to 



subtract from the 'true' map would be a region of immacu- 
late photosphere. Subtract the corresponding level from each 
of the 'difference' or 'scatter' maps. 

(v) Rank each pixel in the 'true' map computed in step 
(iv) against the corresponding pixel in the 'scatter' maps. 
This will show the fraction of fluctuation values that the 
observed value in the 'true' map exceeds. 



For our first test we took the average intensity of the 
outer (and presumably non-irradiated) hemisphere of the 
star as our comparison level, thereby testing the reality of 
the asymmetry between the outer and inner hemispheres. 
The left-hand panel in Fig.|5]shows the resultant significance 
map, calculated using the procedure listed above. The bright 
greyscale depicts regions where the 'true' map ranks above 
100 per cent of the 'scatter' maps, and therefore represents 
regions where the intensity is significantly different from the 
outer hemisphere. Dark greyscales represent regions of the 
map that do not rank highly against the scatter map. Fea- 
tures in these regions are not, therefore, significantly differ- 
ent from the average intensity of the outer hemisphere. 
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Figure 3. Entropy landscapes for the four CVs. Clockwise from top left: AM Her, QQ Vul, HU Aqr and IP Peg. Dark regions indicate 
masses for which no acceptable fit could be found. The cross marks the point of maximum entropy in each landscape. 



The left-hand panel in Fig.|5]shows that the irradiation 
seen on the inner hemisphere in the Roche tomogram is real, 
and appears more significant on the trailing hemisphere (i.e., 
4> = 0.25). Note that the significance of the regions near 
the south pole (SP) is an artefact. Due to the inclination 
of AM Her, latitudes > i in the southern hemisphere are 
never seen and hence the reconstructions always assume the 
(same) default level, with very little scatter in the intensity 
values. As a result, regions near the south pole will appear 
significantly different from the comparison level unless, of 
course, the comparison level is very similar to the default 
level. 

The right-hand panel in Fig. |S] shows another signif- 
icance map, but taking the comparison level as a region 
around the L\ point on the leading hemisphere (<j> — 0.25). 
This shows that the asymmetry between the leading and 
trailing hemispheres is real, and that the effects of irradi- 
ation are primarily found on the trailing hemisphere. The 



asymmetry in the irradiation pattern between the northern 
and southern hemisphere is again an artefact, since regions 
on the southern hemisphere are less visible and the data does 
not constrain the reconstructions in these regions. Note that 
in this case, the comparison level selected was similar to the 
default level of the map, and therefore regions near the south 
pole (SP) no longer appear to be significant. 



Our results are in agreement with the surface maps of 
AM H er obtained using the metho d of radial-velocity curve 
fitting JPavev fc Smithll992l.ri996l) . They proposed that the 
lack of irradiation on the leading hemisphere could be ex- 
plained if the incoming accretion column was able to block 
the radiation produced close to the w hite dwarf and found 
that the accretion geometry given by ICropperl jl98ctl sup- 
ported this argument. 
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AM Her 

primary mass = 0.360 M 
secondary mass = 0.180 M 
orbital period = 11139.00 s 
inclination = 60.00° 
number of slices = 20 
imb darkening OFF 
exposure time = 300.0 s 
intrinsic profile = USER DEFINED 
systemic velocity = 7.00 kms 
instrumental res. = 26.60 kms 
default map = AVERAGE MAP 



X- = 0.70 
entropy = -0.021764 




ocity (km/s) 



velocity (km/s) 



-500 500 

velocity (km/s) 



Figure 4. Roche tomogram of AM Her. Bright greyscales in the tomograms depict regions of enhanced Na I absorption. The grcyscalcs 
are reversed in the trailed spectra. Note that the systemic velocity quoted in the upper-left of the panel does not take into account the 
systemic velocity of the spectral-type template star. 




Figure 5. Significance maps of AM Her. The four images on the left-hand side of each panel are viewed in the equatorial plane, and the 
two images on the right-hand side of each panel are viewed from above and below the north and south pole respectively The significance 
map in the left-hand panel was constructed by taking the comparison level to be the average intensity of the outer hemisphere. The 
significance map in the right-hand panel was constructed by taking the comparison level to be the average intensity of the region around 
the L\ point that does not appear to be irradiated. Bright greyscales indicate regions where the 'true' map ranks above 100 per cent of 
the 'scatter' maps. The significance of the region near the south pole (SP) in the left-hand panel is an artefact, since this region is never 
observed. 
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5 RESULTS FOR QQ VUL 

5.1 The binary geometry of QQ Vul 

The systemic velocity and inclination of QQ Vul were deter- 
mined in the manner outlined in Section 13. 21 and the results 
are shown in Fig.0 We find that the systemic velocity is well 
constrained with a value of -14 km s~ , comparable with 
the va l ue of 17 ± 23 km s -1 obtained bv lMukai fc Charlesl 
Jl987h . lCatalan et al.l (Il999l) found 7 = -0.2 ± 2.2 km s _1 
from the results of elliptical fits to the Na I absorption lines 
used in this work. In addition, using the same method they 
also obtained 7 = -8.3 ± 1.6 km s _1 using data obtained on 
the WHT in 1993. One explanation for this difference in the 
systemic velocity is if a change in the intensity distribution 
between the two epochs can alter the measured value of 7 
(we explore this possibility in Section |7J. 

It is clear from Fig.^that it is difficult to select a single 
value for the inclination given the scatter in the distribution 
of the points. The quality of the reconstructions does, how- 
ever, begin to deteriorate towards lower inclinations allowing 
us to set a lower of limit of ~72° for the inclination of QQ 
Vul. 

Previous estimates of the orbital inclination include i = 



1998 

< 



= 60° ± 14° jMukai & Charles! 11987 


) and 46° < 


jNousek et all |l984|). In addition. 


Nousek et alJ 



( 1984) also found evidence of a brief eclipse of the cy- 
clotron region combined with the simultaneous disappear- 
ance of emission lines. They attributed this to an eclipse 
by the secondary star, forcing the orbital inclination to the 
upper limit of 46° < i < 74°. From EUVE observations 
iBelle. Howell fc Mills! i2000T) also suggested that the X-ray 
emitting region is eclipsed by the secondary star, constrain- 
ing the inclination to > 60° . Certainly, in this work we can 
rule out inclinations < 60° due to the low entropy of the 
reconstructions at these inclinations. These claims of eclipse 
features, however, are outside phase and thus cannot be at- 
tributed to the secondary star. The lack of an X-ray eclipse 
at this phase means that an inclination > 72° can be ex- 
cluded. 

Fig. |3] shows the entropy landscape for QQ Vul, con- 
structed assuming an inclination of 72°, from which we de- 
rive optimum masses of Mi = 0.66 Mq and M2 = 0.42 Mq . 
We adopt these masses and inclination for the reconstruc- 
tion in Section 15.21 The derived masses, however, decrease 
smoothly as the inclination is increased to 90° and cover the 
range Mi = 0.58-0.66 M Q , M 2 = 0.34-0.44 M Q at a constant 
q — 0.63. This agrees well with, and is better constrained 
than, the masses of Mi = 0.54^1 M w and M 2 = 0.30±0.10 
Mq estimated bv lCatalan et alJ |1999). 

The white-dwarf mass w e derive agrees with th e val- 
ues of Mi = 0.5 8±° no Mm JMukai fc Charlesl Il987l) and 



Mi = 0.59 M B (iMouchetl 11993). but is alm ost half that 
of lCropper et alJ lll99Sft and tWu et alJ lll995f) who obtained 
Mi = 1.22 M and Mi = 1.1-1.3 M , respectively. The two 
latter mass determinations make use of the X-rays emitted 
from the accreting white dwarf and involve fitting model 
spectra to the observed X-ray spectra, and are therefore sub- 
ject to underlying assumptions in the model. 

The mass ratio of q = 0.63 obtained in this work is in ex- 
cellent agree ment with previous d eterminations, including q 
= 0.54±0.14 JCatalan et alJll999l) and q = 0.45-0.67 deter- 



mined bv lSchwope et ail (l200Cl) from Doppler tomography 
and polarimetry. 



5.2 The surface map of QQ Vul 

The Roche tomogram of QQ Vul is presented in Fig. |S|and 
shows a distinct reduction in the Na I absorption around 
the L\ point which can be explained by irradiation from 
the white dwarf and accretion regions ionising the Na I on 
the inner face of the secondary, as in AM Her. The effects of 
irradiation also appear to be slightly stronger on the trailing 
hemisphere, suggesting that the leading hemisphere may be 
partially shielded by the accretion stream or curtain. 

Another notable feature is a bright patch on the leading 
hemisphere, corresponding to an increased Na I flux deficit 
around phase 0.75. A similar feature can also be seen in the 
tomogram of IP Peg (Fig. [8J and larger flux deficits at 4> = 
0.75 than at 4> — 0.25 have also bee n reported for the Na 
I doublet in HT Cas JCatalanl Il999l) and in the TiO light 
curves of Z Cha dWade fc HorneHl988T) . 

In their study of QQ Vul u sing the same data used 
in this work, ICatalan et al.l (Il999l) found that their surface 
maps derived from radial-velocity curve fitting did not give 
satisfactory fits to the velocities and line fluxes between or- 
bital phases (j> = 0.6 and <f> — 0.8. The bright patch seen 
around phase 0.75 in Fig. El is almost certainly the cause 
of their unsatisfactory fits. ICatalan et al.l lll999T) attributed 
this discrepancy to the presence of starspots. This interpre- 
tation is almost certainly wrong, however, as an increase 
in the Na I flux deficit is inconsistent with a star-spot, as 
the Na I flux deficit is known t o decrease with later spec- 
tral type jBrett fc Smith! I1993D and star-spots are gener- 
ally cooler than the surrounding photosphere. Therefore, a 
starspot should appear as a dark feature in the tomograms. 
(Although s pots hotter than the photosphere have also been 
miag ed, e.g. lDonati et alJll992l) . 

The reality of the spot feature can be assessed in Fig. |7| 
which shows the significance map when we take the compar- 
ison level as the average intensity on the outer hemisphere, 
avoiding both the region that appears to be irradiated and 
the bright patch. This shows that the bright patch is not 
significant. The significance map does, however, show that 
the irradiated inner hemisphere is real and there is evidence 
for stronger irradiation of the trailing hemisphere. 

As with AM Her, the reduction in absorption around 
the Li point can be explained by irradiation from the white 
dwarf and accreting regions ionising the Na I on the inner 
face of the secondary, with possible shielding by the accre- 
tion stream/curtain reducing the effect on the leading hemi- 
sphere. 



6 RESULTS FOR IP PEG 

6.1 The binary geometry of IP Peg 

The optimum systemic velocity, as shown in Fig. is - 
9 km s _1 . This is significan tly different from the va l ues o f 
56 ± 4 km s" 1 obtained bv iMartin. Jones fc Smith! dl987ft 
and 31 ± 7 km s" 1 obtained bv lMartin et alJ (Il989ft . Fig.[3 
also shows that we have been unable to determine the incli- 
nation of IP Peg. This is due to the small variations in the 
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Figure 6. Roche tomogram of QQ Vul. Bright greyscales in the tomograms depict regions of enhanced Na I absorption. The greyscales 
are reversed in the trailed spectra. Note that the systemic velocity quoted in the upper-left of the panel does not take into account the 
systemic velocity of the spectral-type template star. 



line strength over the orbital cycle, which provides our only 
constraint on the inclination before taking into consideration 
the eclipse width. The reconstructions do, however, yield a 
consistent result of 0.43 for the mass ratio of IP Peg (Fig.|2J. 
This is in disagreement with q = 0.55-0.63 and o = 0.55-0.62 
found bv lMartin et alJ (119891) andlMarsbl <198sh . but agrees 



Crawford 1986) 



0.39T0.04 



with 0.35<q<0 .49 llWood 
jCatalar]ll999T) and a = 0.32±0.08 feeekman et alll200Cl) . 
All the authors, except iBeekman et alJ (120001) . made correc- 
tions for the effects of irradiation where appropriate. 

Although we cannot determine the inclination using the 
surface maps alone, IP Peg is an eclipsing system, and its ge- 
ometry is further constrained by the observed eclipse width, 
which is dependent upon the inclination and size of the sec- 
ondary star's Roche lobe. Fig.[3]shows the entropy landscape 
constructed by varying the inclination at each different mass 
pairing in order to match the eclipse width of A(f> = 0.0863 
iWood fc Crawfordlll986l) . This results in optimum masses 
of Mi = 1.18 M and M 2 = 0.50 M at an inclination of 
82.3°. Altering the eclips e width to the other e xtreme of A(f> 
= 0.0918 given by IWood fc Crawford! dl98rj) changes this 
determination to Mi = 1.16 M Q and M 2 = 0.50 M and 
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Figure 7. The significance map of QQ Vul, taking the compar- 
ison level to be the average intensity of the outer-hemisphere, 
avoiding both the irradiated regions and the bright area on the 
leading hemisphere. Bright greyscales indicate regions where the 
'true' map ranks above 100 per cent of the 'scatter' maps. The 
significance of the region near the south pole (SP) is an artefact 
since this region is never observed. 



the optimum inclination to 84.4°. Our best estimates are 
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Mi = 1.05+° o7 M jBeekman et alJ l2000f) . but slightly 
larger than Mi = 0.85±0.09 M determined by ICatalanl 



- 1 ® 
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Marshl 



and 



trailing hemisphere, but still allow the bright spot to irradi- 
ate the leading hemisphere. 

One further explanation for the features seen around 
the Li point is that they could be caused by the effects 
of an eclipse of the secondary star by the accretion disc. 
This would cause darker regions to appear in the Roche 
tomogram around the hi point and a reduction in the qual- 
ity of the fit. Reconstructions carried out neglecting phases 
between 0.45-0.55 show no improvement in the quality of 
the fit, and there is negligible difference in the Roche to- 
mograms. We are therefore confident that these features are 
not the result of an eclipse. 



6.2 The surface map of IP Peg 

The Roche tomogram of IP Peg in the light of the Na I 
absorption doublet is shown in Fig. |H| The most noticeable 
feature is the weakened absorption flux deficit on the in- 
ner hemisphere of IP Peg, which can be attributed to ir- 
radiation. The effects of irradiation appear to be strongest, 
however, towards the leading hemisphere of the secondary 
star. We also see a bright patch of enhanced Na I absorption 
on the leading hemisphere similar to that seen in QQ Vul 
(Section 

The significance map constructed taking the compari- 
son level to be the average intensity of the outer hemisphere, 
avoiding the bright patch, is shown in Fig. |U] From this we 
can see that the bright patch is not significant, whereas the 
irradiated region on the leading hemisphere is. 

In their radial velocity curve-fitting studies, 
iDavev fc Smith! Jl992T) also found that the effects of 
irradiation were strongest on the leading hemisphere, and 
that regions of IP Peg appeared to be irradiated some 
distance onto the back half of the star. They found that 
this feature could not be explained by irradiation from 
the bright spot, as the illuminating source would have 
to originate from a direction about 45° away from the 
line of centres and inst ead attributed this as ymmetry to 
circulation currents (e.g. IMartin fc Davevl ll995'). 

It is clear in Fig. [S] that we cannot identify any signifi- 
cant irradiation effects on the back half of the star and, as 
such, cannot provide any support for the evidence of circu- 
lation currents on IP Peg. It seems likely, instead, that the 
irradiation of the leading hemisphere can be explained by 
the bright spot, which is located on the correct side of the 
secondary star to cause this effect, an d is a strong illumi- 
nating source (as seen in, for example. ISzkodvlll987l) . 

Our Roche tomogra m, however, bears l ittle r esemblance 
to the surface map that IDavev fc Smith! dl992l) computed 
from a model for a secondary star heated equally by the 
bright spot and the white dwarf (their Fig. 15), in that there 
does not appear to be any significant signs of irradiation on 
the trailing hemisphere. This suggests that irradiation from 
both the white dwarf and the bright spot is shielded from 
this region. An obvious way to block irradiation from the 
surface of the white dwarf is through shielding by the accre- 
tion disc. Irradiation from the bright spot may be blocked 
by the accretion stream, especially if much of the emission 
from the bright spot is located downstream from the im- 
pact point, as s ugges ted in the eclipse maps of IP Peg by 
iBobinger et alJ (Il999l) . This would result in shielding of the 



7 RESULTS FOR HU AQR 

7.1 The binary geometry of HU Aqr 

Fig. shows the systemic velocity and inclination of HU 
Aqr as a function of entropy. We obtain a value of 7 = 
7 km s _1 and, although it is difficult to assign a value to 
the inclination, the reconstruction quality begins to decrease 
rapidly for inclinations less than 80° , and we therefore adopt 
a best estimate of i > 80°. 

Our estimate of the systemic velocity is c ompatible at 
the 2g level with 7 = 83 ± 44 km s _1 derived bv lGlenn et alJ 
lll994h who fitted a sinusoid to the narrow Ha component 
thoug ht to originate from the secondary star. lSchwope et alJ 
( 1997) derived a systemic velocity of -7.8 ±1.4 km s _1 by 
fitting radial velocity curves to the same data used in this 
work. Although this result disagrees with the value obtained 
using Roche tomography, the sine-fitting method takes no 
account of any non-uniform intensity distribution present on 
the secondary star, whereas Roche tomography does. Since 
there is clearly an asymmetrical intensity distribution on the 
secondary star in HU Aqr (see Section f7.2H we believe that 
this is the cause of the difference between the two results. 

In order to test the reliability of the sine-fitting method 
in determining the systemic velocity, we generated a model 
secondary star with the same binary geometry as we 
have found for HU Aqr in this work (Mi = 0.61 M , 
M 2 = 0.15 M and i = 84.4°, see later) but with 7 = 
km s _1 . The inner hemisphere of the model was then uni- 
formly irradiated except for a totally dark shadow covering 
40 per cent of the leading hemisphere. The outer hemisphere 
of the secondary star was also set to be totally dark. From 
this model, trailed spectra of infinite signal-to-noise were 
generated, assuming 100 equally spaced phase bins and 30 
second exposure lengths. 

By fitting radial-velocity curves to this fake data, we 
obtained a systemic velocity of -14 km s _1 . Fig. 1101 shows 
the result of the radial-velocity curve fit (solid line) to 
the synthetic trailed spectrum, and the true motion of 
the centre-of-mass of the secondary star (dashed line) for 
comparison. The value of 7 = -14 km s _1 is significantly 
different from the true value of 7 = km s _1 used in 
the model, and brings our systemic velocity of 7 km s _1 
into agre ement with the value of -7.8 ± 1.4 km s _1 ob- 
tained bv lSchwope et al] lll997l) . This shows that a knowl- 
edge of the intensity distribution across the secondary is 
required to ac curately determine t he systemic velocity (as 
also found by ISchwope et alJ I1997T) . This may have impli- 
cations for the accurate determinations of 7's used in CV 
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Figure 8. Roche tomogram of IP Peg. Bright greyscales in the tomograms depict regions of enhanced Na I absorption. The greyscalcs 
are reversed in the trailed spectra. Note that the systemic velocity quoted in the upper-left of the panel does not take into account the 
systemic velocity of the spectral-type template star. 




= 0.25 




(f> = 0.5 



= 0.75 




Figure 9. The significance map of IP Peg, taking the comparison 
level to be the average intensity of the outer-hemisphere avoiding 
the bright region on the leading hemisphere. Bright greyscalcs 
indicate regions where the 'true' map ranks above 100 per cent 
of the 'scatter' maps. 



Paradiis. Aueusteiin fc S tchlc 199^ and 



age tests fe.g..|van 
iNorth et al.ll2002h . 

Since HU Aqr is an eclipsing system, its inclination has 
been reasonably constrained by previous studies, including 



estimates of i = 80° ±5 ° jGlenn et al.lll994Tl obtained from 
fits t o polarisation curves, and i = 85.6° iSchwope et alJ 
l200ll) from modelling the soft X-ray light-curves. Our deter- 
mination of i > 80° does not provide any further constraint 
to those values already published. 

Previous estimates of the mass ratio, however, are 
somewhat more uncertain and r ange from roughly 0.18 
- 0.4 (e.g. ISchwope et alJ Il99l lHakala et all Il993l and 
iHeerlein. Home fc Schwopd ll999'l, Fig.[5]shows the compo- 
nent masses and mass ratio obtained when reconstructions 
are carried out over different inclinations. This yields a mass 
ratio for HU Aqr of q = 0.25, w hich is in excellen t agreement 
with q = 0.256 derived bv IHeerlein et af assuming 1 

= 84°, and the value of q = 0.25 that lSchwope et alJ (1200 ll) 
found. 

Fig. |3 shows the entropy landscape for HU Aqr which 
has been constructed by changing the inclination at each 
different mass pa iring in order to ma tch the eclipse width 
of A<f> = 0.0763 iHeerlein et al.lll999l) . as described for IP 
Peg in Section 16.11 This results in a best estimate of the 
component masses of Mi = 0. 61 My and Mi = 0.1 5 Mq. 
This agrees with the result of ISchwope et ail J2001), who 
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Figure 10. The synthetic trailed spectrum used to test the reli- 
ability of the sine-fitting method for the purposes of determining 
the systemic velocity. See the text for a description of the model 
used to generate the synthetic data ("Section 17.11 . The solid line 
corresponds to the radial-velocity curve fit and the dashed line 
corresponds to the true motion of the centre-of-mass of the sec- 
ondary star. Note that the radial-velocity curve fit has also ad- 
justed the zero point to <j> = 0.02. 

determined M2 = 0.15 Mousing the mass-radius relation of 
ICaillault fc Patterson! lll990l). a nd M 2 = 0.17 M according 
to the relation bv lNeeca <ll984f) . Our derived masses, com- 
bined with the eclipse width, give our best estimate of the 
inclination as 84.4°. 



7.2 The surface map of HU Aqr 

The Roche tomogram of HU Aqr in the light of the He II 
4686A emission line is presented in Fig. 1111 The two most 
prominent features in the Roche tomogram are the strong 
asymmetry between the inner and outer hemispheres of the 
secondary star and the weaker asymmetry between the lead- 
ing and trailing hemispheres. The emission from the inner 
hemisphere is due to irradiation from the primary and the 
accretion regions, and the tomogram is compatible with 90 
per cent of the emission arising from the inner hemisphere 
alone. (Note that the remainder of the emission that arises 
on the outer hemisphere is most likely due to the smearing 
effects of the default map). 

The reality of the asymmetry between the leading and 
trailing hemispheres can be assessed from the left-most panel 
in Fig. 1121 which shows a slice passing through the leading 
hemisphere (LH), north pole (NP), trailing hemisphere (TH) 
and south pole (SP) near the L\ point of the secondary star. 
The triangular points represent the 'true' map values; 67 
per cent of the bootstrapped 'trial' maps (measured relative 
to the mode of the distribution) lie within the solid line, 
and 100 per cent within the dotted line 1 . It is clear from 
this that the asymmetry between the leading and trailing 
hemispheres is not due to statistical errors. This asymmetry 

1 We have not generated a significance map for HU Aqr since the 
low scatter in the pixel values derived from the bootstrap recon- 
structions mean that each region of the map would be significant 
when compared with any other region. 



was also found bv lSchwope et all <ll997f) who attributed it 
to shielding by the accretion curtain. 

Closer examination of the tomogram reveals more sub- 
tle intensity variations across the secondary star. In partic- 
ular there appears to be a region of more intense emission 
located towards the north pole but shifted slightly towards 
the leading hemisphere. The plot on the right hand-side of 
Fig. H2l shows a vertical slice cutting through the secondary 
star close to the polar regions but slightly towards the L\ 
point. This again shows the shadow on the leading hemi- 
sphere (LH) as discussed earlier, but another feature of note 
is the drop in emission from the trailing hemisphere (TH) 
compared to the polar regions. This feature is not apparent 
near the L\ point and only becomes significant for regions 
halfway between the centre-of-mass of the secondary and 
the L\ point, disappearing again when the effects of irra- 
diation are diminished near the terminator. This cannot be 
explained by shielding of the secondary star by the accre- 
tion curtain, since this would only block irradiation on the 
leading hemisphere. 

HU Aqr does, however, have a strong stream component 
with as much as 40 per cent of the flux in the U band due 
to accretion stream emission jHakala et al] ll993L Doppler 
maps bv lSchwope et alJ Jl997h also show that the accretion 
stream is a dominant emission-line feature. Presumably the 
accretion stream can also irradiate the secondary star and, 
since the stream is deflected in the direction of the orbital 
motion, the stream will be invisible to regions of the sec- 
ondary star on the trailing hemisphere away from the L\ 
point. This obscuration of the accretion stream by the sec- 
ondary star itself may explain the observed shadow on the 
trailing hemisphere. 

This feature could also be explained if systematic errors 
were able to produce a greater line flux in the centres of the 
lines, corresponding to artificially brighter polar regions in 
the Roche tomograms and giving the illusion of a shadow at 
lower latitudes. It was first thought that the assumption of 
a Gaussian profile for the emission from the secondary star 
could introduce such an artefact, since the actual profile 
from the secondary star is expected to be flatter, especially 
around quadrature. We found, however, that at low spectral 
resolutions the emission profile from the secondary is well 
approximated by a Gaussian. Artefacts could be generated 
if the emission from the various components of HU Aqr (e.g. 
from the secondary star and the accretion stream) are in- 
correctly disentangled, which is especially difficult at phases 
where the individual components are crossing each other. 
This situation is more difficult to assess and emphasises the 
importance of mapping lines that are known to originate 
solely from the secondary star. 



8 CONCLUSIONS 

The Roche tomograms show an asymmetric irradiation pat- 
tern on all 4 CVs. This can be explained if the accretion 
column/stream shields the leading hemisphere in AM Her, 
QQ Vul and HU Aqr. There is also evidence for a fainter 
shadow present on the trailing hemisphere of HU Aqr. We 
propose that the stream itself contributes significantly to the 
irradiation, and that the shadow on the trailing hemisphere 
occurs as these regions on the secondary star have no direct 
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Figure 11. Roche tomogram of HU Aqr. Bright grcyscalcs in the tomograms depict regions of enhanced He II A4686A emission. The 
greyscales are reversed in the trailed spectra. 




100% 

95% 

— 67% " 



95% 
67% 




Figure 12. The left-hand panel shows a vertical slice (sec text for explanation) through HU Aqr near the L\ point; the right-hand panel 
shows a vertical slice through HU Aqr closer to the terminator. The triangular points show the intensities found in the original fit. The 
lines show confidence limits derived from 200 bootstrapped reconstructions. The labels on the top axis denote the positions of pixels 
located on the leading hemisphere (LH), north pole (NP), trailing hemisphere (TH) and south pole (SP). 
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view of the majority of the stream. Further observations us- 
ing absorption lines originating from the secondary star are 
needed to confirm this. 

A different scenario occurs in IP Peg, where irradiation 
affects the leading hemisphere only. We suggest that the 
bright spot, which is on the correct side of the secondary star 
to cause such an irradiation pattern, may be the principal 
illuminating source in this system. Although we still expect 
irradiation of the trailing hemisphere by the white dwarf (as 
seen in the magnetic CVs), there is no evidence for this in 
the tomogram, ft may be the case that the accretion disc in 
IP Peg is shielding the secondary star f r om ir radiation by 
the white dwarf. Unlike ll)avev fe Smith! il992T) . we find no 
evidence for irradiation extending onto the back half of the 
star, and hence cannot support their idea that circulation 
currents operate in IP Peg. 

We have used the entropy landscape technique to de- 
rive the system parameters. Although we cannot reliably 
constrain the inclination, and therefore the masses, in the 
non-eclipsing CVs, the technique does provide a consistent 
result for the mass ratio that is independent of the incli- 
nation. For the two eclipsing CVs (IP Peg and HU Aqr), 
we believe that the masses derived in this work are the 
most accurate to date. In addition, it is found that measure- 
ments of the systemic velocity using circular orbit fits to the 
radial-velocity variations are unreliable, and are dependent 
upon the intensity distribution across the secondary star. We 
might therefore expect such measurements to be both line 
and time-dependent. Since the intensity distribution across 
the secondary star is known in Roche tomography, systemic 
velocity measurements using this technique should prove far 
more accurate. 
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